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SUMMARY 

The abrupt transition to heightened poliomyelitis epidemicity in England and Wales, 1947-1957, 
was associated with a profound change in the spatial dynamics of the disease. Drawing on the 
complete record of poliomyelitis notifications in England and Wales, we use a robust method of 
spatial epidemiological analysis (swash-backwash model) to evaluate the geographical rate of 
disease propagation in successive poliomyelitis seasons, 1940-1964. Comparisons with earlier and 
later time periods show that the period of heightened poliomyelitis epidemicity corresponded with 
a sudden and pronounced increase in the spatial rate of disease propagation. This change was 
observed for both urban and rural areas and points to an abrupt enhancement in the propensity 
for the geographical spread of polioviruses. Competing theories of the epidemic emergence of 
poliomyelitis in England and Wales should be assessed in the light of this evidence. 
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The same general sequence of poliomyelitis emergence 
(sporadic cases small outbreaks—*- major epidemic 
events) was repeated during the late nineteenth and 
early twentieth centuries in other countries of 
Europe, North America and elsewhere [1, 2, 10]. By 
the 1950s, poliomyelitis had assumed the status of 
an epidemic disease of global proportions [1, 2, 11]. 

The emergence of epidemic poliomyelitis in 
twentieth-century England and Wales assumed a dis- 
tinctive form, with a step-function transition to heigh- 
tened epidemicity in the early post-war years (Fig. la). 
Between the world wars, the monthly rate of polio- 
myelitis notifications in England and Wales rarely 
rose above 0-5 cases/100000 population; 1926 and 
1938 were the only years of epidemic note [12, 13]. 
The situation changed suddenly and markedly there - 

* Author for correspondence: Dr M. Smallman-Raynor, School of after. Beginning with the first nationwide epidemic 
Geography, University of Nottingham, University Park, . n « r . , juri , j j ji 

Nottingham NG7 2RD, UK. m 1947 ' En S land and Wales entered a decade-long 

(Email: matthew.smallman-raynor@nottingham.ac.uk) phase of heightened poliomyelitis epidemicity [2]. 

The online version of this article is published within an Open Access environment subject to the conditions of the Creative Commons 
Attribution-NonCommercial-ShareAlike licence <http://creativecommons.Org/licenses/by-nc-sa/3.0/>. The written permission of Cambridge 
University Press must be obtained for commercial re-use. 



INTRODUCTION 

Epidemic poliomyelitis was one of the important 
emergent viral diseases of the twentieth century 
[1, 2], While several lines of evidence point to the 
occurrence of poliomyelitis in pre-modern societies 
[3-5], medical references to the disease prior to the 
nineteenth century are few in number and are con- 
sistent with the circulation of an endemic infection 
associated with sporadic cases of clinical disease [1]. 
Reports of small and highly localized outbreaks 
of 'infantile paralysis' began to gain momentum 
in Europe from the 1880s. The world's first major 
epidemics occurred in Scandinavia in 1905 [6-9]. 
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Fig. 1. Monthly series of poliomyelitis notification rates per 100000 population, 1913-1971. (a) England and Wales 
(source: General Register Office, London); (b) USA (source: Public Health Reports/ Morbidity and Mortality Weekly 
Report). 



This emergence profile differs from the experience of 
the USA where, in the wake of the epidemic of 
1916, a relatively smooth and progressive transition 
to the major epidemic events of the post-war years 
was apparent (Fig. lb) [14, 15]. On both sides of the 
Atlantic, the post-war period of heightened epidemi- 
city was brought to a rapid close by the introduction 
and widespread administration of safe and effective 
poliovirus vaccines from the mid-1950s [1, 2]. 

The epidemic emergence of poliomyelitis in post- 
war England and Wales was accompanied by a pro- 
nounced change in the geography of the disease 
[14], with the replacement of small focal outbreaks 



(inter-war years) by large national epidemics (post- 
war years) signalling a fundamental shift in the spatial 
dynamics of disease propagation [12, 16]. An in-depth 
analysis of the geographical spread of the 1947 
epidemic appears elsewhere [17]. In the present article, 
we examine the shifting pattern of poliomyelitis 
activity with particular reference to the geographical 
rate of disease propagation in the 25-year sequence 
of poliomyelitis seasons, 1940-1964. Our examination 
draws on a robust method of spatial epidemiological 
analysis known as the swash-backwash model [18-20]. 
This model is a spatial derivative of the generic 
SIR (susceptible =>■ infected =>■ recovered) mass action 
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Table 1. Mean and range of mid-point population estimates for the 1479 standard local government districts of 
England and Wales at the beginning and end of the 25-year observation period, 1940-1964* 







Mean (range)t 




Category of district 


n 


1940 


1964 


Urban areas 








Metropolitan boroughs 


28 


109-8 (23-5-271-3) 


113-6 (20-4-348-5) 


County boroughs 


83 


152-2 (23-3-997-9) 


164-0 (31-6-1102-8) 


Municipal boroughs 


310 


28-3 (0-9-156-5) 


33-7 (1-0-182-6) 


Urban districts 


576 


13-8 (0-8-188-7) 


181 (0-5-210-3) 


Rural areas 








Rural districts 


482 


16-0 (0-4-78-6) 


20-8 (1-5-104-2) 



* Data from Registrar General's Statistical Review of England and Wales (London: HMSO). 
t Population estimates in thousands. 



models commonly used to model the spread of infec- 
tious diseases in human and many animal populations 
[21, 22]. Using the binary presence/absence of a dis- 
ease rather than the actual number of notified cases 
in a set of districts, the swash-backwash model: 
(i) allows the identification of the leading and trailing 
edges of the spatial advance ('swash') and retreat 
('backwash') of an infection wave; (ii) provides a 
means of measuring the phase transitions of geo- 
graphical units from susceptible (S), through infective 
(I) and recovered (R) status as dimensionless integrals; 
and (iii) yields robust measures of the spatial velocity 
of epidemic waves from space-time series. 

Our analysis demonstrates that the period of 
heightened poliomyelitis epidemicity in England and 
Wales, 1947-1957, was associated with a sudden and 
pronounced increase in the geographical rate of dis- 
ease propagation. The change occurred in all types 
of geographical areas, affecting large, medium and 
small urban centres and rural districts alike. 



METHODS 
Data 

Statutory case notifications of poliomyelitis in 
England and Wales have been collated by the 
General Register Office, London, since 1 September 
1912 [23]. For the purposes of the present analysis, 
we drew on the weekly notifications of poliomyelitis 
(acute poliomyelitis and acute polioencephalitis) 
included in the Registrar General's Weekly Return of 
Births and Deaths . . . and Cases of Certain Specified 
Infectious Diseases (London: HMSO) for 1940-1964. 
For each of the 1305 registration weeks (ending 



on a Saturday) in the 25-year observation period, 
information on the number of notified cases of 
poliomyelitis in the 1479 standard local government 
districts (28 metropolitan boroughs, 83 county 
boroughs, 310 municipal boroughs, 576 urban dis- 
tricts and 482 rural districts; see Table 1) of England 
and Wales was abstracted from the Weekly Return 
to form a (1305 weeks x 1479 districts=) 1-93 million 
cell matrix of poliomyelitis notifications. This time- 
space matrix included records for a total of 69928 
cases of poliomyelitis. 

The matrix cells for the constituent districts of each 
of the (then) 62 counties of England and Wales were 
grouped to form composite 1305-week time series of 
poliomyelitis notifications for each of three geographi- 
cal divisions of a given county: (i) metropolitan and 
county boroughs, representing large urban areas; 
(ii) municipal boroughs and urban districts, represent- 
ing medium and small urban areas; and (iii) rural 
districts (Table 1). Some categories of district were 
absent from some counties, so that the aggregation 
procedure resulted in a total of 152 county-specific 
time series: metropolitan and county boroughs, 
31 counties; municipal boroughs and urban districts, 
61 counties; and rural districts, 60 counties. The set 
of 1 52 time series (representing an exhaustive division 
of England and Wales) and the subsets of 31, 61 and 
60 time series (representing different levels in the 
settlement system) form the basis of our analysis. 

Defining the poliomyelitis season, 1940-1964 

The seasonal proclivity of poliomyelitis is one of the 
outstanding epidemiological features of the disease 
[24-26]. During the inter- and early post-war years, 
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Fig. 2. The poliomyelitis season, England and Wales, 1940-1964. For a given year (1940,..., 1964), the count of 
poliomyelitis notifications in each registration week, w (1,..., 53) was assessed relative to the median weekly count 
of notifications for that year; weeks above the median were defined as 'excess weeks'. The bar chart plots the number of 
years in which registration week w was categorized as an 'excess week' in the 25-year observation period. The median 
value of the set of bars (8) is demarcated by the broken horizontal line; registration weeks above this median are defined 
as constituting the poliomyelitis season (registration weeks 24-50). The heavy line trace is shown for reference and plots 
the median poliomyelitis notification rate per 100000 population for registration week it' in the period 1940-1964. 



the months of June-December formed the typical 
season of raised poliomyelitis activity in the temperate 
latitudes of the northern hemisphere [2]. To define 
formally the poliomyelitis season for England and 
Wales, 1940-1964, we adopted a robust method that 
apportioned equal weight to the temporal signal in 
each of the 25 observation years. The method, which 
is summarized in Figure 2, defines the poliomyelitis 
season as disease registration weeks 24-50 (mid-June 
to mid-December). Weekly time series of poliomyelitis 
notifications are plotted for successive poliomyelitis 
seasons in Figure 3; the aggregate count of notifi- 
cations in each season, along with information on 
the predominant circulating virus type(s) isolated 
from patients with clinical disease, appears in Table 2. 



The swash-backwash model 

To assess the geographical rate of poliomyelitis 
propagation in each of the poliomyelitis seasons, 
1940-1964, we used the swash-backwash model of 
spatial epidemiological analysis [18-20]. In essence, 
the model measures the rate of spatial advance and 
retreat of an infection wave across a set of geographi- 
cal units (e.g. the counties of England and Wales or 
the states of the USA). Let the first week of a given 
poliomyelitis season (registration week 24) be coded 
as t - 1 . The subsequent weeks of each season are 



then coded serially as t = 2, 3, . . ., 27. For each county, 
we refer to the first week of the season in which a 
case of poliomyelitis was notified in any one of the 
three categories of districts (metropolitan and county 
boroughs; municipal boroughs and urban districts; 
and rural districts) as the leading edge (LE) and the 
last week in which a case was notified as the following 
edge (FE) of the season in that category. A series of 
measures of the spatial velocity of an infection wave 
can be derived from the LE and FE and these are 
defined in the Technical Appendix. Summary descrip- 
tors of the measures in relation to the present analysis 
are given in Table 3. 

Average velocity measures (Vle, Vee, Ife — ^le) 

For each edge, LE and FE, we can define a time- 
weighted arithmetic mean, F LE and f FE , which gives 
the average time in weeks of (notified) arrival (LE) 
and departure (FE) of the infection wave in a 
given category of districts across the set of counties 
[Technical Appendix, equation (Al)]. These time- 
weighted means can be converted to dimensionless 
velocity ratios F LE and V FE with values in the range 
[0, 1] [equation (A2)]. In addition, the difference 
between F F e and Fle (i-e. Ffe - Fle) provides a measure 
of the average duration of the poliomyelitis wave in a 
given category of district (Table 3). We refer to this 
measure as the average duration of infectivity. 
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Fig. 3. Time series of poliomyelitis notifications, England and Wales, 1940-1964. Graphs plot the weekly count of notifications 
for successive poliomyelitis seasons (registration weeks 24-50). 



Integrals (S A , I A , R A ) 

As described in the Technical Appendix, the integrals 
S A [the proportion of the study area A at risk of infec- 
tion; equation (A4)], I A [the proportion of A which is 
infected; equation (A5)], and R A [the proportion of A 
which has recovered; equation (A6)] can be derived 
from the average measures of spatial velocity / LE and 
Ffe- AH three integrals are dimensionless numbers with 
values in the range [0, 1] such that S A + I A + R A -\. 
As summarized in Table 3, low values of S A denote 
a very rapid spatial advance of an infection wave 
while high values of R A denote very rapid spatial 
retreat. Finally, I A provides a measure of the relative 
rates of spatial advance and retreat of an infection 
wave. 

Spatial basic reproduction rate ( Rq A ) 

In conventional mass action (SIR) epidemic models, 
the basic reproduction rate, Rq, is one of the most 



useful parameters for the mathematical characteriz- 
ation of infectious disease processes [equation (A3)] 
[22, 29]. The spatial analogue of R 0 , the spatial basic 
reproduction rate R 0A , is derived from the integrals 
S A and R A [equation (A8)] and provides a measure 
of the propensity of an infected geographical unit to 
spawn other infected units in later time periods. 
Effectively, values of Rq A calibrate the spatial velocity 
of disease spread, with higher values denoting more 
rapid spatial propagation (Table 3). 

Model application 

Based on the epidemiological curve in Figure la, the 
25-year observation interval (1940-1964) was sec- 
tioned into three major time periods: (i) 1940-1946, 
representing the interval prior to the period of heigh- 
tened epidemicity; (ii) 1947-1957, bracketed by the 
first (1947) and last (1957) major epidemic years and 
representing the period of heightened epidemicity; 
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Table 2. Application of the swash-backwash model to the 1940-1964 poliomyelitis seasons, England and Wales* 





Predominant 


Poliomyelitis 














Average duration 


Seasonj 


poliovirus type| 


notifications 


S A 


1a 


Ha 






RoA 


of infectivity§ 


1940 


? 




889 


0-38 


0-39 


0-23 


0-59 


0-23 


2-69 


9-60 


1941 


? 




699 


0-37 


0-38 


0-24 


0-59 


0-24 


2-56 


9-28 


1942 


? 




518 


0-36 


0-36 


0-28 


0-61 


0-28 


2-30 


8-83 


1943 


? 




323 


0-40 


0-27 


0-33 


0-56 


0-33 


1-82 


6-30 


1944 


? 




405 


0-31 


0-34 


0-35 


0-65 


0-35 


1-97 


8-13 


1945 


? 




699 


0-39 


0-41 


0-21 


0-58 


0-21 


2-97 


1001 


1946 






552 


0-32 


0-41 


0-27 


0-64 


0-27 


2-49 


9-99 








OOJZ 


A.I o 


U-OV 


ft. i "> 


ft.T7 
U - / / 


U-1Z 


(TOO 


1 /*o4 


!V4o 




i 


IjjD 


a.*>s 

U Zo 


U JJ 


ft.l o 
U IV 


ft./^O 


ft. 10 
U If 


J ol 


1 j jy 


1 0 i 0 








u zu 




ft.i i 

U 11 


ft.TT 
U* / / 


ft.l 1 
U 11 


T.ft^ 

/ Uj 


1 T.^ft 


1950 






7753 


0-15 


0-74 


0-11 


0-81 


0-11 


7-61 


18-89 


1951 


i, m 


2453 


0-17 


0-62 


0-21 


0-79 


0-21 


3-87 


15-63 


1952 


i, m 


3788 


0-18 


0-65 


016 


0-78 


016 


4-97 


16-58 


1953 


i,n 


4392 


014 


0-68 


018 


0-82 


018 


4-85 


17-47 


1954 


i 




1820 


0-22 


0-59 


0-20 


0-75 


0-20 


4-00 


14-87 


1955 


i, m 


6755 


0-22 


0-67 


011 


0-75 


011 


6-95 


17-09 


1956 


i 




3104 


019 


0-61 


0-20 


0-77 


0-20 


4-13 


15-53 


1957 


i 




4272 


016 


0-63 


0-21 


0-80 


0-21 


4-08 


1613 


1958 


i 




1828 


0-25 


0-53 


0-22 


0-71 


0-22 


3-45 


13-29 


1959 


i 




879 


0-25 


0-42 


0-33 


0-72 


0-33 


2-30 


10-47 


1960 


i, in 


368 


0-31 


0-33 


0-36 


0-65 


0-36 


1-94 


8-03 


1961 


i, in 


854 


0-40 


0-35 


0-26 


0-57 


0-26 


2-36 


8-40 


1962 


i, in 


220 


0-24 


0-23 


0-53 


0-72 


0-53 


1-43 


5-14 


1963 






51 


0-37 


018 


0-45 


0-59 


0-45 


1-39 


3-80 


1964 






28 


0-53 


011 


0-36 


0-44 


0-36 


1-31 


2-00 



* The period of heightened poliomyelitis activity, 1947-1957, is shown in bold type, 
t Registration weeks 24-50. 

| Based on information included in annual editions of the Report of the Ministry of Health (London: HMSO), with sup- 
plementary data from Goffe [27] and Spicer [28]. 
§ Formed as Fee — 'le (measured in weeks). 



and (iii) 1958-1964, representing a period of epidemic 
retreat associated with the introduction and mass 
administration of safe and effective poliovirus vac- 
cines. The swash-backwash model was then applied 
to the sets of poliomyelitis seasons (registration 
weeks 24-50) in each of the three time periods. 
Analysis was undertaken for the entire set of 1 52 geo- 
graphical divisions and for each of the three subsets 
of divisions (metropolitan and county boroughs; 
municipal boroughs and urban districts; and rural 
districts). All model parameters were calculated on 
a cumulative (weekly) basis for each poliomyelitis 
season, thereby permitting an assessment of velocity 
as the season unfolded. 

Data issues and model interpretation 

The principal issues encountered in the geographical 
analysis of notified cases of poliovirus infection 



have been reviewed elsewhere [2]. Here, we note 
that the characteristically high ratio of subclinical 
to clinical cases of poliovirus infection [30] dictates 
that only a fraction of the total number of infec- 
tions in the population of England and Wales will 
have been captured in the Registrar General's 
Weekly Return over the 25-year observation period, 
1940-64. Under these circumstances, the swash- 
backwash model provides a measure of the spatial 
rate of advance and retreat of notified cases of clinical 
poliovirus infection. The use of notified cases of 
infection as a proxy for the underpinning infection 
wave is subject to a number of potentially confound- 
ing factors, including spatial and temporal variations 
in the subclinicalxlinical infection ratio and the 
completeness of the reporting of clinically apparent 
infections [2, pp. 53-59]. Our interpretation of 
the results which follows is subject to these data 
caveats. 
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Table 3. Summary descriptors of the parameters of the swash-backwash model as applied to poliomyelitis seasons in 
England and Wales, 1940-1964 

Model 

parameter Description 

K LE Dimensionless velocity ratio in the range [0, 1], measuring the average time from the onset of the poliomyelitis 

season to the first notified case of poliomyelitis in a given category of district. Relatively low (high) values denote 

relatively slow (fast) spreading infection waves. 
Vpe Dimensionless velocity ratio in the range [0, 1], measuring the average time from the onset of the poliomyelitis 

season to the last notified case of poliomyelitis in a given category of district. Relatively low (high) values denote 

relatively slow (fast) retreating infection waves. 
'fe — 'le Average duration of infectivity, defined as the difference (in weeks) between the average time to the first (/le) and 

last (Ffe) notified cases of poliomyelitis in a given category of district. Relatively low (high) values denote 

infection waves of relatively short (long) duration. 
S A Susceptible integral, defines the proportion of geographical units in the study area A which are at risk of infection 

with poliomyelitis. Relatively low (high) values denote relatively fast (slow) spreading infection waves. 
I A Infected integral, defines the proportion of geographical units in the study area A which are infected with 

poliomyelitis. The integral forms a measure of the relative rate of spatial advance and retreat of an infection wave. 

Relatively low (high) values denote infection waves in which the rate of spatial retreat is relatively fast (slow) 

compared with the rate of spatial advance. 
R\ Recovered integral, defines the proportion of geographical units in the study area A which have recovered from 

infection with poliomyelitis. Relatively low (high) values denote relatively slow (fast) retreating infection waves. 
Roa Spatial basic reproduction rate, measures the propensity of an infected geographical unit to spawn other infected 

units in later time periods. Values of Rq A calibrate the spatial velocity of poliomyelitis spread, with relatively low 

(high) values denoting relatively slow (fast) spatial propagation. 



RESULTS 

Aggregate analysis: all geographical divisions 

Figure 4 plots the cumulative percentage of the 
152 geographical divisions that first notified (leading 
edge, LE) and last notified (following edge, FE) 
cases of poliomyelitis in each registration week of 
the poliomyelitis season. Cumulative percentages 
were created as mean estimates across the correspond- 
ing registration weeks for 1940-1946, 1947-1957 and 
1958-1964. The shaded envelopes define the 95% 
confidence intervals for the estimates. The graphs 
depict the mean (spatial) phase transitions of the 152 
geographical divisions from susceptible to infected 
(S =^ I) (Fig. 4a) and infected to recovered (I =>■ R) 
(Fig. 4b) status in the three time intervals. 

Taken relative to the averages for earlier 
(1940-1946) and later (1958-1964) poliomyelitis 
seasons, Figure 4 highlights two key features of 
the spatial transmission of poliomyelitis in the 
period of heightened epidemicity (1947-1957): (i) the 
much more rapid growth of the LE of the infection 
wave (Fig. 4a), indicative of the relatively faster 
spatial advance of poliomyelitis; and (ii) the slower 
growth of the FE of the infection wave (Fig. 4b), 
indicative of the relatively slower spatial retreat of 
poliomyelitis. 



Swash-backwash parameters 

Figures 5 and 6 plot, in the manner of Figure 4, 
weekly mean estimates of the integrals S A , I A and 
R A (Fig. 5a-c), the dimensionless velocity parameters 
Fle and F FE (Fig. 6a, b), the spatial basic repro- 
duction rate, R 0A (Fig. 6c) and the average duration 
of infectivity ?fe — ?le (Fig. 6d) for the poliomyelitis 
seasons of 1940-1946, 1947-1957 and 1958-1964. 
As before, the shaded envelopes define the 95% confi- 
dence intervals for the mean estimates. Parameter 
values, calculated to the end of the poliomyelitis 
season (registration week 50), are given for individual 
years in Table 2 and as mean estimates (with 95% 
confidence intervals) in the upper rows of Table 4. 

The substantially lower values of S A and higher 
values of F LE confirm the more rapid spatial advance 
of poliomyelitis in 1947-1957 than in the preceding 
(1940-1946) and following (1958-1964) periods 
(Tables 2 and 4). Similarly, the substantially lower 
values of R A and V FE confirm the slower spatial 
retreat of poliomyelitis in 1947-1957 than in the 
earlier and later time periods (Tables 2 and 4). 
Figures 5 and 6 show that inter-period differences in 
parameter values become evident at approximately 
3-10 weeks from the beginning of the average polio- 
myelitis season, with 1947-1957 being differentiated 
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Fig. 4. Spatial phase transitions of poliomyelitis. The solid lines plot, as mean values, the cumulative proportion of 
districts that (a) first notified (leading edge, LE) and (b) last notified (following edge, FE) cases of poliomyelitis in each 
week of the poliomyelitis seasons (registration weeks 24-50) for 1940-1946, 1947-1957 and 1958-1964. The shaded 
envelopes define the 95% confidence levels for the mean estimates. The graphs compare the evidence for 1947-1957 with 
1940-1946 (left) and 1958-1964 (right). 



from the other time periods by the maintenance of a 
considerably faster rate of spatial advance (S A , F LE ; 
Figs 5a, 6a) and a considerably slower rate of spatial 
retreat (R A , F FE ; Figs 5c, 6b) from thereon. The inter- 
period differences are further reflected in the much 
higher values of I A for 1947-1957 (Tables 2 and 4, 
Fig. 5b), indicative of a relatively more rapid phase 
transition from susceptible to infected (S =>■ I) than 
infected to recovered (I =>■ R) status in the years of 
heightened epidemicity. 

Overall, the higher rate of spatial propagation of 
poliomyelitis in 1947-1957 is underscored by the 
higher values of R 0A (Tables 2 and 4, Fig. 6c), indi- 
cative of the greater propensity for infected geographi- 
cal units to yield secondaries in this time period. 
Inspection of Figure 6c shows that differences in the 
values of R 0A were most pronounced in the early 
and intermediate stages of an average poliomyelitis 
season, with the values for the three time periods 
converging towards the end of the season. 

The average duration of infectivity (Ffe — ?le) high- 
lights inter-period differences in the average length of 



time that cases of poliomyelitis were reported from 
geographical units (Fig. 6d), with much more durable 
infection waves in 1947-1957 (>16 weeks) than 
1940-1946 and 1958-1964 (<9 weeks) (Table 4). 

Finally, viewed on a year-by-year basis, Table 2 
shows that the onset of heightened epidemicity 
(1947) was associated with an abrupt change in the 
values of all the model parameters, signifying a sud- 
den shift in the geographical rate of disease propa- 
gation that was maintained for the next 10 years. 

Disaggregated analysis: urban and rural areas 

Figure 7 and Table 4 summarize the results of the 
swash-backwash analysis for each of the three subsets 
of geographical divisions: metropolitan and county 
boroughs; municipal boroughs and urban districts; 
and rural districts. For each subset of geographical 
units, Figure 7 plots the average weekly values of 
the LE (graph a), FE (graph b) and the swash- 
backwash parameters (graphs c-i) for the poliomyel- 
itis seasons of 1940-1946, 1947-1957 and 1958-1964. 
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Fig. 5. Swash-backwash parameters, I. The solid lines plot, for each week of the poliomyelitis season (registration weeks 
24-50), mean values of (a) S A , (b) I A and (c) R A for 1940-1946, 1947-1957 and 1958-1964. The shaded envelopes define 
the 95% confidence levels for the mean estimates. The graphs compare the evidence for 1947-1957 with 1940-1946 (left) 
and 1958-1964 (right). 



Average values of the swash-backwash parameters, 
formed to the last week of the poliomyelitis season 
(registration week 50), are given for each subset of 
divisions in the lower rows of Table 4. The findings 
of the aggregate analysis are repeated for each 
subset of geographical divisions in Figure 7 and 
Table 4, indicating that the change in spatial 
dynamics associated with the period of heightened 
poliomyelitis epidemicity (1947-1957) affected large 
urban areas, medium and small urban areas and 
rural areas alike. 



DISCUSSION 

Theories of the mechanisms that may have promoted 
the epidemic emergence of poliomyelitis in the 
twentieth century are numerous and embrace a 
range of biological, environmental, social and demo- 
graphic considerations [1, 2, 31], While changes in 
the human host, the viral agent and the virus trans- 
mission mechanism have all been invoked in an under- 
standing of the emergence complex [2], the role of 
hygiene as a factor in shifting patterns of population 
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Fig. 6. Swash-backwash parameters, II. The solid lines plot, for each week of the poliomyelitis season (registration weeks 
24-50), mean values of (a) V LE , (b) V FE , (c) R 0A and (d) duration of infectivity (? F e-?le) for 1940-1946, 1947-1957 and 
1958-1964. The shaded envelopes define the 95% confidence levels for the mean estimates. The graphs compare the 
evidence for 1947-1957 with 1940-1946 (left) and 1958-1964 (right). 



exposure and immunity has received particular atten- 
tion in the literature [24, 31-33]. According to the 
hygiene model, successive improvements in levels of 
hygiene and sanitation during the late nineteenth 
and early twentieth centuries account for a reduced 
level of faecal exposure to poliovirus in early infancy, 
thereby reducing the level of latent immunization in 



the population. Further perspectives on the model 
are provided elsewhere [24, 33] but, underpinned by 
progressive developments in sanitation, the hygiene 
model would account for the gradual development 
of more sizable epidemics of poliomyelitis, associated 
with older patient cohorts, as the twentieth century 
progressed. 
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Table 4. Mean values of the swash-backwash parameters for the 1940-1964 poliomyelitis seasons* 

















Average duration 


Period 


Sa 


Ia 


Ra 


He 


v FE 


Rqa 


of infectivityt 



All geographical divisions (n = 152) 

1940-1946 0-36 (±0-03) 0-37 (±0-04) 0-27 (±0-04) 0-60 (±0-02) 0-27 (±0-04) 2-40 (±0-30) 8-88 (±0-98) 

1947-1957 0-19 (±0-02) 0-65 (±0-03) 0-16 (±0-03) 0-77 (±0-02) 0-16 (±0-03) 5-29 (±0-89) 16-44 (±0-91) 

1958-1964 0-33 (±0-08) 0-31 (±0-11) 0-36 (±O08) 0-63 (±0-08) 0-36 (±0-08) 2-02 (±0-56) 7-30 (±2-91) 

Metropolitan and county boroughs (n = 31) 

1940-1946 0-31 (±0-04) 0-44 (±0-08) 0-25 (±0-05) 0-65 (±0-04) 0-25 (±0-05) 2-93 (±0-63) 10-82 (±2-15) 

1947-1957 0-14 (±0-02) 0-75 (±0-02) 0-11 (±0-02) 0-83 (±0-02) 0-11 (±0-02) 8-49 (±1-61) 19-26 (±0-62) 

1958-1964 0-28 (±0-11) 0-43 (±0-15) 0-29 (±0-10) 0-68 (±0-11) 0-29 (±0-10) 3-03 (±1-22) 10-56 (±4-12) 

Municipal boroughs and urban districts (« = 61) 

1940-1946 0-36 (±0-02) 0-38 (±0-03) 0-27 (±0-04) 0-61 (±0-02) 0-27 (±0-04) 2-47 (±0-25) 9-18 (±0-72) 

1947-1957 0-19 (±0-03) 0-65 (±0-03) 0-16 (±0-02) 0-77 (±0-03) 0-16 (±0-02) 5-23 (±0-62) 16-47 (±0-84) 

1958-1964 0-34 (±0-09) 0-32 (±0-11) 0-34 (±0-08) 0-63 (±0-09) 0-34 (±0-08) 2-10 (±0-55) 7-72 (±2-85) 

Rural districts (n = 60) 

1940-1946 0-40 (±0-03) 0-32 (±0-04) 0-29 (±0-05) 0-57 (±0-03) 0-29 (±0-05) 2-19 (±0-36) 7-53 (±1-06) 

1947-1957 0-21 (±0-02) 0-60 (±0-05) 0-19 (±0-04) 0-75 (±0-02) 0-19 (±0-04) 4-65 (±1-09) 15-13 (±1-26) 

1958-1964 0-36 (±0-05) 0-21 (±0-10) 0-43 (±0-10) 0-60 (±0-05) 0-43 (±0-10) 1-62 (±0-44) 4-67 (±2-77) 

* Parameter values formed to the final week of the poliomyelitis season (week 50). 95% confidence intervals for the mean esti- 
mates are given in parentheses, 
t Formed as f FE — f LE (measured in weeks). 



The specific factors that promoted the emergence of 
epidemic poliomyelitis in England and Wales, and 
which could account for the sustained period of heigh- 
tened epidemicity with abrupt onset in 1947, have yet 
to be fully deciphered. Some contemporary observers 
highlighted the perceptible shift in the incidence of 
poliomyelitis to older age groups which, on the basis 
of fragmentary information, had begun to manifest 
in England and Wales in the mid- 1930s and which is 
consistent with the hygiene model of poliomyelitis 
emergence [12, 14, 34, 35]. Available evidence, how- 
ever, indicates that the increase in the age of attack 
was not gradual but was characterized by a series 
of abrupt shifts (notably, in the 'epidemic' year of 
1938 and, again, in the post-war years) and this 
has prompted speculation over the possible role of 
additional emergence mechanisms [12]. Agerholm, 
for example, attributed the step-function increase 
in poliomyelitis activity (Fig. la) to the introduction 
of some 'new' poliovirus type, possibly Brunhilde 
(type 1), for which the British population had limited 
immunity [36]. There is insufficient information avail- 
able on the poliovirus types that were circulating in 
Britain prior to the 1950s to test this hypothesis, 
although we note that poliovirus type 1 was the pre- 
dominant virus type identified in a small sample of 
clinical cases that occurred in the period 1946-1950 
(Table 2) [27]. 



Framed by these perspectives, our analysis of polio- 
myelitis data for England and Wales has shown that, 
relative to the immediately preceding (1940-1946) and 
following (1958-1964) years, the period of heightened 
poliomyelitis epidemicity (1947-1957) was associated 
with: (i) a faster rate of spatial advance of polio- 
myelitis; (ii) a slower rate of spatial retreat of polio- 
myelitis; and (iii) an extended period of notified 
poliomyelitis activity. These changes were under- 
pinned by a shift in the geographical pattern of disease 
activity, from small focal outbreaks (inter-war years) 
to national epidemics (post-war years) [14], and are 
consistent with the operation of an emergence process 
that considerably enhanced the efficiency by which 
poliomyelitis spread from one geographical area to 
another. Importantly, this enhancement of geographi- 
cal spread was not gradual; it occurred abruptly with 
the onset of heightened epidemicity and involved both 
urban and rural areas of the country. 

Cognisant of the apparent limitations of the 
hygiene model as the sole explanatory device for the 
transition to heightened poliomyelitis epidemicity in 
post-war England and Wales [17], our analysis points 
to the role of additional factors that may have con- 
tributed to the enhanced geographical efficiency of 
disease transmission at the time. In particular, we 
note that the diffusion characteristics of the events 
of 1947-1957 are analogous to those observed when 
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Fig. 7. Swash-backwash parameters for urban and rural areas. The graphs plot, for each week of the poliomyelitis season 
(registration weeks 24-50), mean values of the swash-backwash parameters for 1940-1946, 1947-1957 and 1958-1964. 
The information is shown for metropolitan and county boroughs (left), municipal boroughs and urban districts (centre) 
and rural districts (right). 



the swash-backwash model is applied to the first-time observation is consistent with the hypothesized intro- 
spread of viral diseases such as new (pandemic) strains duction of a virulent poliovirus type into England 
of influenza A in a population [20]. While this and Wales in the early post-war period [36], other 
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possible explanations (including the prevalence of 
unusually transmissible poliovirus strains in the 
years from 1947) merit careful consideration. 

While the documented increase in poliomyelitis 
activity in post-war England and Wales was a real 
epidemiological phenomenon and not an artefact of 
disease reporting [37-39], caution should be exercised 
in the interpretation of the poliomyelitis-related infor- 
mation included in the Registrar General's Weekly 
Return. First, as we have already noted, the high 
ratio of subclinical to clinical cases of poliovirus 
infection [30] dictates that only a fraction of the 
total number of infections in the population of 
England and Wales will have been captured by the 
notification system. A corollary of this data limitation 
is that the observed geographical patterns of polio- 
myelitis activity, described in this paper, may differ 
in unknown ways from the underpinning patterns of 
poliovirus infection. Second, although assessments of 
the inter-area comparability of poliomyelitis notifi- 
cations for England and Wales at mid-century are 
relatively favourable [34], unknown spatial and tem- 
poral biases in the notification of severe (paralytic) 
and less severe (non-paralytic) disease are likely to 
exist. Third, all disease notifications included in the 
Weekly Return are provisional and are subject to 
correction upon receipt of additional information 
from local health officers. Corrected data, however, 
are only available on a temporally aggregated 
(quarterly and annual) basis and we have followed 
the Ministry of Health in our use of (uncorrected) 
weekly notifications in the present article [38]. The 
swash-backwash model has provided us with a robust 
method for examining the spatial dynamics of polio- 
myelitis under these conditions of uncertain data 
quality. The more general utility of the model for 
the robust spatial analysis of infectious disease events 
is emphasized. 



departure of the poliomyelitis wave in the set of 
geographical units. For the LE, the equation is 

1 T 

?LE = ^X/"" (A1) 



(=1 



where n t is the number of units whose LE occurred 
in epidemic week t and N—^n t . The time-weighted 
mean is a useful measure of the velocity of the wave 
in terms of average time to unit infection. A similar 
equation can be written for FE, and higher-order 
moments can also be specified. To allow comparison 
between geographical areas with epidemic waves 
of different duration, these can be converted into a 
velocity ratio, V (0s= Fs= 1) 



Vi 



LE 



1 - 



^LE 



(A2) 



where T is the duration of the wave. A similar 
equation can be written for K FE . 

Spatial basic reproduction rate (R 0A ) 

In conventional SIR epidemic models, the basic repro- 
duction rate, R 0 , is defined as the ratio between an 
infection rate (/?) and a recovery rate (y): 

Ro=P/y. (A3) 

In the spatial domain, A, the spatial basic reproduction 
rate, R 0A , is the average number of secondary infected 
geographical units produced from one infected unit in 
a virgin area. In a given study area, the integral S A 
(the proportion of the study area at risk of infection) 
is given by 

(Zle-1) 
T 



(A4) 



while the proportion of the area which is infected 
(the infected area integral) is 



?FE 



-s A . 



(A5) 



TECHNICAL APPENDIX 

Following Cliff & Haggett [18], let the weeks of an 
epidemic be serially coded t=\, 2, 3,..., T. For a 
given geographical unit, we refer to the first week in 
which a case of acute poliomyelitis is notified as the 
leading edge (LE) and the last week in which a case 
is notified as the following edge (FE) of the infection 
wave in that district. For each edge we can define 
a time-weighted arithmetic mean, ? LE and ?fe This 
gives us the average time of (notified) arrival and 



The recovered area integral, R A , is 
R A = l-(S A + I A ). 



(A6) 



Rearranging equation (A5) and substituting into 
equation (A6), R A can be expressed as 



Ra = 1 



(A7) 



Thus, R A is equivalent to the dimensionless velocity 
measure V FE . 

The integrals in equations (A4)-(A6) are dimen- 
sionless numbers with values in the range [0, 1]. 
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The integral S A has parallels to the reciprocal of fl in 
that a small value of S A suggests very rapid spread. 
The integral R A has parallels to y in that a large 
value suggests very rapid recovery. Since S A is inver- 
sely related to /?, Cliff and colleagues [18, 40] proposed 
a spatial version of R 0 which, for the purposes of the 
present analysis, takes the form: 

R 0A = 1 -^. (A8) 

As Cliff & Smallman-Raynor [40, pp. 147-178] 
observe, other formulations of R 0A are possible and 
have been used with equal success [18-20]. The spatial 
reproduction number in equation (A8) measures the 
propensity of an infected geographical unit to spawn 
other infected units in later time periods. In effect, 
values of R 0A calibrate the velocity of such spread 
(the larger the value, the greater the rate of spread). 

It is important not to over-stretch the analogy 
between R 0 and R 0A . Ro is defined for the hypothetical 
situation when a new case is introduced into a wholly 
susceptible population. While R 0A is defined for a 
virgin region, it is calculated using spatial data for 
the entire span of the outbreak. As a result, it is 
contaminated with data from the later phases of the 
outbreak when many spatial sub-areas are no longer 
virgin. This may account for the frequent small calcu- 
lated values. Ro is useful because it distinguishes 
between situations where an epidemic can take off 
(Ro>l) and those where it cannot (Rq<1), and this 
is arguably the most important attribute of R 0 as a 
summary parameter in epidemiology. Rq's spatial 
cousin, R 0A , does not share this property - Cliff 
et al. [20] provide examples of real epidemics that 
had sustained spread from one district to another 
in which R 0A < 1 . But it does allow large numbers 
of spatial settings to be examined and the relative 
velocity at different stages of outbreaks to be assessed 
and compared. Finally, it should be noted that the 
parameters R 0A , LE and FE are correlated, but each 
gives slightly different insights into the progress of 
an outbreak through a geographical area. 

ACKNOWLEDGEMENTS 

We are grateful to the anonymous referees for their in- 
sightful comments on an earlier version of this paper. 

DECLARATION OF INTEREST 

None. 



REFERENCES 

1. Paul JR. The History of Poliomyelitis. New Haven: 
Yale University Press, 1971, pp. 486. 

2. Smallman-Raynor M, et al. Poliomyelitis. A World 
Geography: Emergence to Eradication. Oxford: Oxford 
University Press, 2006, pp. 734. 

3. Hamburger O. A case of infantile paralysis in antiquity 
[in French]. Bulletin de la Societe Francaise d'Histoire de 
la Medecine 1911; 10: 407^112. 

4. Henschen F. The History of Diseases. London: 
Longman, 1966, pp. 344. 

5. Dequeker J, Fabry G, Vanopdenbosch L. Hieronymous 
Bosch (1450-1516): paleopathology of the medieval dis- 
abled and its relation to the Bone and Joint Decade, 
2000-2010. Israel Medical Association Journal 2001; 3: 
864-871. 

6. Wickman I. Beitrage zur kenntnis der Heine-Medinschen 
krankheit (poliomyelitis acuta und verwandter erkran- 
kungen). Berlin: S. Karger, 1907, pp. 292. 

7. Holt LE, Bartlett FH. The epidemiology of acute polio- 
myelitis. A study of thirty-five epidemics. American 
Journal of the Medical Sciences 1908; 135: 647-662. 

8. Batten FE. The epidemiology of poliomyelitis. Proceed- 
ings of the Royal Society of Medicine 1911; 4: 198-226. 

9. Low RB. A short account of the epidemiology of acute 
anterior poliomyelitis in recent years. Forty-fifth Annual 
Report of the Local Government Board, 1915-16. Supple- 
ment Containing the Report of the Medical Officer for 
1915-16. London: HMSO, 1917, pp. 14-74. 

10. May JM. Map of the world distribution of poliomyeli- 
tis. Geographical Review 1950; 40: 646-648. 

1 1 . Nathanson N, Kew OM. From emergence to eradica- 
tion: the epidemiology of poliomyelitis deconstructed. 
American Journal of Epidemiology 2010; 172: 1213-1229. 

12. Bradford Hill A. Poliomyelitis in England and Wales 
between the wars. Proceedings of the Royal Society of 
Medicine 1954; 47: 795-806. 

13. MacNalty A. Poliomyelitis. In: Cruickshank R, ed. 
Control of the Common Fevers. London: The Lancet 
Ltd, 1942, pp. 292-303. 

14. Martin WJ. Poliomyelitis in England and Wales, 
1947-50. British Journal of Social Medicine 1951; 5: 
236-246. 

15. Trevelyan B, Smallman-Raynor M, Cliff AD. The 

spatial dynamics of poliomyelitis in the USA: from epi- 
demic emergence to vaccine-induced retreat, 1910-1971. 
Annals of the Association of American Geographers 
2005; 95: 269-293. 

16. Hardy A. Poliomyelitis and the neurologists: the view 
from England, 1896-1966. Bulletin of the History of 
Medicine 1997; 71: 249-272. 

17. Smallman-Raynor MR, Cliff AD. The geographical 
spread of the 1947 poliomyelitis epidemic in England 
and Wales: spatial wave propagation of an enigmatic 
epidemiological event. Journal of Historical Geography 
2013; 40: 36-51. 

18. Cliff AD, Haggett P. A swash-backwash model of the 
single epidemic wave. Journal of Geographical Systems 
2006; 8: 227-252. 



Poliomyelitis in England and Wales, 1947-1957 591 



19. Cliff AD, Haggett P, Smallman-Raynor MR. An 

exploratory method for estimating the changing speed 
of epidemic waves from historical data. International 
Journal of Epidemiology 2008; 37: 106-112. 

20. Cliff AD, et al. Infectious Diseases. Emergence and Re- 
emergence: A Geographical Analysis. Oxford: Oxford 
University Press, 2009, pp. 763. 

21. Bailey NT J. The Mathematical Theory of Infectious 
Diseases and Its Applications. London: Griffin, 1975, 
pp. 413. 

22. Anderson RM, May RM. Infectious Diseases of 
Humans: Dynamics and Control. Oxford: Oxford 
University Press, 1992, pp. 766. 

23. Local Government Board. Forty-second Annual Report 
of the Local Government Board, 1912-13. Supplement 
Containing the Report of the Medical Officer for 
1912-13. London: HMSO, 1914. 

24. Melnick JL. Poliovirus and other enteroviruses. In: 
Evans AS, Kaslow RA, eds. Viral Infections of 
Humans: Epidemiology and Control, 4th edn. London: 
Plenum Medical Book Company, 1997, pp. 583-663. 

25. Dowell SF. Seasonal variation in host susceptibility and 
cycles of certain infectious diseases. Emerging Infectious 
Diseases 2001; 7: 369-374. 

26. Fisman DN. Seasonality of infectious diseases. Annual 
Review of Public Health 2007; 28: 127-143. 

27. Goffe AP. Studies on the distribution of poliomyelitis 
viruses in England and Wales. Proceedings of the 
Royal Society of Medicine 1955; 48: 937-941. 

28. Spicer CC. The incidence of poliomyelitis virus in 
normal children aged 0-5 years. A report on a 
study by the Public Health Laboratory Service and 
Local Health Authorities. Journal of Hygiene 
(London) 1961; 59: 143-159. 

29. Heffernan JM, Smith RJ, Wahl LM. Perspectives on the 
basic reproductive ratio. Journal of the Royal Society 
Interface 2005; 2: 281-293. 

30. Heymann DL. Control of Communicable Diseases 
Manual, 18th edn. Washington, D.C.: American 



Public Health Association, 2004, pp. 700. 

31. Paul JR. Epidemiology of poliomyelitis. In: Polio- 
myelitis. Geneva: World Health Organization, 1955, 
pp. 9-29. 

32. Sabin AB. Epidemiologic patterns of poliomyelitis in 
different parts of the world. In: Poliomyelitis: Papers 
and Discussions Presented at the First International 
Poliomyelitis Conference. Compiled and Edited for the 
International Poliomyelitis Congress. Philadelphia: Lippin- 
cott, 1949, pp. 3-33. 

33. Nathanson N, Martin JR. The epidemiology of 
poliomyelitis: enigmas surrounding its appearance, 
epidemicity and disappearance. American Journal of 
Epidemiology 1979; 110: 672-691. 

34. Daley A. Acute anterior poliomyelitis in 1947 with 
special reference to London. Proceedings of the Royal 
Society of Medicine 1948; 41: 52-60. 

35. Benjamin B, Gale AH. The age distribution of polio- 
myelitis in England and Wales. Monthly Bulletin of 
the Ministry of Health and the Public Health Laboratory 
Service 1949; 8: 208-213. 

36. Agerholm M. Epidemic poliomyelitis and its prevention. 
Lancet 1953; 265: 287-293. 

37. Murray LH. The 1947 outbreak of poliomyelitis: an 
interim report to the week ended 20th September, 
1947. Monthly Bulletin of the Ministry of Health and 
the Public Health Laboratory Service 1947; 6: 169-173. 

38. Gale AH. Poliomyelitis in England and Wales in 
1947. Monthly Bulletin of the Ministry of Health and 
the Public Health Laboratory Service 1948; 7: 127-132. 

39. Ministry of Health. Report of the Ministry of Health for 
the Year Ended 31st March 1948, Including the Report of 
the Chief Medical Officer on the State of the Public 
Health for the Year Ended 31st December 1947. 
London: HMSO, 1949. 

40. Cliff AD, Smallman-Raynor MR. Oxford Textbook of 
Infectious Disease Control: A Geographical Analysis from 
Medieval Quarantine to Global Eradication. Oxford: 
Oxford University Press, 2013, pp. 208. 



